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Much experimental and theoretical efforts have been devoted in the past twenty years to search for 
a genuine thermodynamic reentrant phase transition from a ferromagnetic to either a paramagnetic 
or spin glass phase in disordered ferromagnets. So far, no real system or theoretical model of a 
short-range spin glass system has been shown to convincingly display such a reentrant transition. 
We present here results from Migdal-Kadanoff real-space renormalization-group calculations that 
provide for the first time strong evidence for ferromagnetic to paramagnetic reentrance in Potts 
spin glasses on hierarchical lattices. Our results imply that there is no fundamental reason ruling 
out thermodynamic reentrant phase transitions in all non mean-field randomly frustrated systems, 
and may open the possibility that true reentrance might occur in some yet to be discovered real 
randomly frustrated materials. 



I. INTRODUCTION 

All real materials contain a certain amount of frozen-in 
random disorder. Often, random disorder leads to ran- 
domly competing, or frustrated, interactions jy. Ran- 
dom frustration is detrimental to the type of order 
that an otherwise idealized pure material would dis- 
play for zero disorder. Randomly frustrated systems are 
ubiquoutous in condensed matter physics. Examples in- 
clude: magnetic systems ||, mixed molecular crystals ||, 
superconducting Josephson junctions in an applied mag- 
netic field 0], liquid crystals in porous media [||, and 
partially UV polymerized membranes l|. 

One of the main issues at stake in all frustrated systems 
is how the low temperature phase of the pure material 
is affected by weak disorder, and how that state evolves 
with increasing disorder level. In particular, one of the 
most intriguing questions is whether a weakly frustrated 
system can lose upon cooling the long-range ordered 
phase established at higher temperature and return, or 
reenter, into a thermally disordered phase, or go into a 
randomly frozen glassy phase. Because of their relative 
simplicity over other systems, and because of the large 
number of systems readily available with easily control- 
lable level of disorder, random magnets are ideal systems 
to study the effects of weak frustration and to investigate 
the above question. It was originally thought that sev- 
eral weakly frustrated ferromagnets, such as Eu^Sri-^S 
and amorphous-Fei_ x Mna;, were displaying a reentrant 
transition from a fcrromagnetically long-range ordered 
phase to a randomly frozen spin glass phase upon cool- 
ing and for a finite range of disorder, x Q. However, af- 
ter twenty years of extensive experimental research, it is 
now generally believed that a true thermodynamic reen- 
trant phase transition from a long-range ordered ferro- 
magnetic (F) phase to either a spin glass (SG) or para- 
magnetic (P) phase does not occur in real magnetic ma- 



terials ■ Once established at the P-F Curie tem- 
perature, T c , ferromagnetic order remains down to zero 
temperature, though with the possibility of a transverse 
spin-freezing transition at Tj_ (0 < Tj_ < T c ) in XY 
and Hciscnberg systems, which does not destroy the fer- 
romagnetic order Above a critical disorder level, 
ferromagnetism occurs only on short length scales, and 
the system displays, instead, full-blown spin-glass behav- 
ior below a glass transition temperature, T g 

At the theoretical level, it is also currently believed 
that reentrance does not occur in any random bond spin 
glass model J7|-[2l|. This is certainly the case for the in- 
finite range Ising, XY and Heisenberg models [^. In 
two and three dimensional Ising and Heisenberg models, 
high-temperature series expansion Monte Carlo 

simulations (^0,^1 and recent defect-wall energy calcu- 
lations ||lj find no reentrant behavior either. Recently, 
compelling renormalization-group p~2|Jl~3) and quenched 
gauge symmetry arguments 16 Jl7]havc been put for- 
ward for a broad class of spin glass models, which include 
the Ising spin glass [^) and the gauge glass model for dis- 
ordered Josephson junction arrays and vortex glass in dis- 
ordered type-II superconductors p| Jl9| , |20| , p^ | , and which 
strongly argue against reentrance. Some of the details 
of these predictions have been quantitatively tested by 
high-temperature series expansion (ll] , while possibly 
the most detailed checks have been obtained from real- 
space Migdal-Kadanoff renormalization group (MKRG) 
calculations of Ising spin glass models on so-called hi- 
erarchical lattices [p~0| — p~2| - Even in the case of the 
two-dimensional XY model with random Dzyaloshinskii- 
Moriya couplings, which for a long-time was believed to 
be a good candidate for reentrant behavior (|,|2j|, evi- 
dence is now rapidly accumulati ng th at reentrance does 
not occur in that system either |p[ ^9| , ^o| . 

Summing up, it appears that the case against reen- 
trance in randomly frustrated systems and non mean- 
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field theoretical models p4| is at this time simply over- 
whelming (2^]. In fact, the evidence is sufficiently strong 
that it could be interpreted as an indication that some 
profound, though yet unknown, reason(s) formally forbid 
reentrance in all spin glasses, even those which do not ex- 
hibit a quenched gauge invariance |pl1 , |l2| , ^6|Jl7| l . This is 
not impossible given that our understanding of the nature 
of the ground state(s) and of the low-lying excitations in 
glassy systems is still limited. In this letter we present a 
counter example to this common belief, what we believe 
is the first strong evidence for a reentrant transition in 
a simple non mean-field spin glass model where thermal 
fluctuations and the question of lower-critical spatial di- 
mension play a key role ^4j. Specifically, we consider 
the 3-state ferromagnetic Potts model with a concentra- 
tion x of random antiferromagnetic bonds on hierarchi- 
cal lattices. We investigate the thermodynamic behavior 
of this model using the MKRG scheme, which is an ex- 
act method for hierarchical lattices We present re- 
sults which show that this model exhibit ferromagnetic to 
paramagnetic reentrance for a finite concentration range 
of antiferromagnetic bonds. Reentrance is made possible 
by the fact that the system prefers to lower its free en- 
ergy through short range antiferromagnetic (AF) corre- 
lations rather than to preserve long-range ferromagnetic 
order. Since the lower critical dimension for antiferro- 
magnetic long range order for the Q = 3 Potts model 
on hierarchical lattices is 4 ^tJ, the system can be reen- 
trant in 2 and 3 dimensions. It is interesting to note 
that the MKRG method has in the past been used as 
one of the key methods in establishing the absence of 
reentrance in Ising |^fj|- |l2| , ^4| , XY || and possibly also 
in Dzyaloshinskii-Moriya XY spin glasses |) 19 20|p3] |. 



II. MODEL AND METHOD 



The Hamiltonian for the Q-state Potts model is: 



H 



E * 



(i) 



<i,3> 



where Jij > for ferromagnetic couplings and < 
for antiferromagnetic ones. 5 ai ,a- is the Kronecker delta: 
the spin, <n at lattice site i can take Q states, Q — 
0, 1, 2, . . . Q — 1. The bond energy between two spins 
is — Jij if the two spins are in the same state cr^ = Uj, 
and zero otherwise. The familiar Ising model is equiva- 
lent to a Q = 2 state Potts model with a shift of total 
energy of the system, and a rescaling of the exchange cou- 
pling by a factor 2. Although "less popular" than the 
Ising model, the 3-state Potts model is also important 
in modelling real condensed matter systems. For exam- 
ple, the two-dimensional antiferromagnetic 3-state Potts 
model on the frustrated kagome lattice captures some of 
the essential of the low temperature thermodynamics of 
the Heisenberg antiferromagnet on that lattice p8| . Also, 
it has been suggested that the orientational freezing in 



molecular glasses, such as N2-Ar and KBr-KCN, can be 
partially described by a three-dimensional 3-state Potts 
spin glass model || . 

Here we consider the situation where the bonds Jij in 
Eq. ^ are distributed randomly, and given by a quenched 
biased bimodal probability distribution, V(Jij): 



V{Jij) = xS(J tJ - J) + (1 - x)8{J lJ + J) 



(2) 



A bond between sites i and j has a probability x to be 
ferromagnetic and of strength J, and a probability 1 — x 
of being AF and of strength — J. We study the thermo- 
dynamic properties of this system on hierarchical lattices 
using the MKRG scheme §^|l|Jl|,|l|,||,|o) . One consid- 
ers a sequence of b Jy bonds in series, each we label 
jW (k = 1,2,. ..6), where (6-1) spins are summed over 
(we have dropped the subscript ij). The above Hamil- 
tonian preserves its invariant form (apart from a spin- 
independent term) under the decimation of the (6 — 1) 
spins. This results in a new effective coupling Jij(l + 1), 
at the RG decimation step 
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k=b 
k=l 



{ CX p{/3J(fc>(i)}-l} 
(3) 

and (3 = 1/k^T. In dimension d, M* 4-1 ) such parallel 
paths of b bonds in series, each with its end-to-end cou- 
pling J n (l + 1), are then added together to give one new 
coupling 



Jij(l + 1) 



E 

71=1 



J (n) (^ + 1) 



(4) 



In practice, the procedure is implemented by first cre- 
ating a large pool of N (N w 10 6 ) bare couplings, 
Jij (I — 0), distributed according to Eq. ||. Then, b cou- 
plings are randomly picked out of that pool combined to 
create a serial coupling J n (l + 1) as given byEq. |. Then, 
b^ -1 ) such couplings J n (l + 1) are added together to give 
one new coupling J^ (1+ 1). The procedure is repeated N 
times to repopulate a new pool of N couplings Jij (I + 1) 
at RG step {I + 1). 

The nature of the magnetic phase at a given tempera- 
ture, T, and concentration, x, of antiferromagnetic bonds 
is determined by monitoring the I— dependence of the av- 
erage value, J(l), and the width, AJ(l), of the distribu- 
tion of N bonds Jij (l). As I — > 00, J and AJ evolve in 
the various phases as 



lim J —> , lim A J — > 

I— >oo / — >-oo 



,• 7 r AJ 

urn J — * +00 , lim — =- 

l — >OG I — >OG J 

lim J — > —00 , lim — ^~ 

l — >oo / — >oc J 

lim A J — > 00 , lim — — - 

/ — >oo / — >oo Zaj 



paramagnetic (5) 

: ferromagnetic (6) 

: antiferromagnetic (7) 

: spin glass (8) 
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To allow for the existence of an antiferromagnetic phase, 
we must work with odd values of 6, as even values of 
b "frustrate" the antiferromagnetic phases and maps an 
initial startup antiferromagnetically biased V(Jij) into a 
ferromagnetic phase already at iteration step #1. Here 
we focus on hierarchical lattices with b = 3. 



III. RESULTS 

The temperature vs concentration of antiferromagnetic 
bonds phase diagram for the three-dimensional d = 3 
case (with b — 3) is shown in Fig. 1. The phases have 
been determined according to the criteria given above. 
Firstly, there is no AF or SG phase at nonzero temper- 
ature in this system in the whole range < x < 1. The 
most remarkable feature of this phase diagram is the exis- 
tence of a reentrant ferromagnetic to paramagnetic phase 
transition for the range x £ [0.765,0.855]. The value of 
0.855 obtained by extrapolating these nonzero tempera- 
ture results agrees with the one obtained by iterating the 
MKRG equations above at zero temperature exactly [^9| . 
Similar results were obtained for d = 2. 
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FIG. 1. Temperature, T, concentration, x, of antiferromag- 
netic bonds phase diagram for the Q — 3 state Potts spin glass 
onad = 3,6 = 3 hierarchical lattices. Ferromagnetic to para- 
magnetic reentrance occurs in the range 0.765 < x < 0.855. 
There is no spin glass or antiferromagnetic phase in this model 
at nonzero temperature. 



As observed 
discussed in other papers MM 



and 

|19|,^9|J30| , the MKRG 
scheme is difficult to implement for spin glass models at 
low temperatures (as T — > 0+). The reason for this is as 
follows. The occurence of a ferromagnetic phase is mon- 
itored by the criterion [J — * +oo, AJ/J — > 0]. In prac- 
tice the ferromagnetic phase is detected when numerical 
overflow occurs as I — ► oo on the ferromagnetic side of 
the P— F boundary. In a spin glass model, where the dis- 
tribution of exp{(3Jij(l)} is broad at low-temperature in 
the ferromagnetic phase close to either the P— F or F— SG 



boundary, one encounters overflow at iteration step / max 
way before J(l) has increased by several order of mag- 
nitude compared to J (I = 0). In previous MKRG stud- 
ies of spin glasses |-y|,|Ul|] the F-P or the F-SG 
phase boundary did not give any "peculiar" reentrant 
boundary, and the extrapolated finite temperature F— P 
or F— SG boundary down to T = agreed with explicit 
MKRG calculations at T = 0. Consequently, there has 
been until now no incentive to push the limit of the nu- 
merics in MKRG calculations of spin glasses as T — > 0+. 
However, in our case here, with this novel reentrant be- 
havior, one could be concerned that the lower reentrant 
portion of the phase boundary is a numerical artifact. 
Specifically, it would a priori seem possible that the flow 
between 0.765 and 0.855 seems to indicate a paramag- 
netic phase according to the criterion given above, for 
1 < I < Imax, but actually, be found to "reverse itself" 
for a value I > l TCV with Z rcv > Z max , if numerical overflow 
bounds allowed it to be seen, and such that the asymp- 
totic large length scale behavior for 0.765 < x < 0.855 
was ferromagnetic in the limit I — > oo. In such a sce- 
nario, the reentrant region would result from a combina- 
tion of short length scale physics added to a finite limit 
to overflow bounds imposed by the computer used for the 
calculations. 

To address this issue, we parametrized each of the 
exp{/3Jij(l)} "coupling terms" via a two-component vec- 
tor exp{/3Jij(l)} = {Mijify^ijQ)}, where Mij(l) and 
Eij(l) are the mantissa and the exponent, in base 10, of 
exp{/3 Jij (I) . The MKRG computer code was then rewrit- 
ten in terms of direct algebraic mantissa operations and 
exponent shifting operations. With this improved version 
of the MKRG computer code, the upper limit for over- 
flow for double-precision calculations on a 32-bit machine 
moves from 10 308 to ~ 10^ 10 a tremendious improve- 
ment. With this modification, the MKRG iterations be- 
come for all practical purposes devoid of overflow limita- 
tions. Our results with this version of the MKRG scheme 
gave an identical phase boundary to the one obtained 
using straighforward conventional double-precision cal- 
culations on a 32-bit machine. The results in Fig. 1 for 
T/J< 0.25 were actually obtained with the "improved" 
version of the MKRG scheme. We are therefore confi- 
dent that the reentrant phase transition displayed by the 
Q = 3 bimodal Potts spin glass model on the b = 3 hi- 
erarchical lattice is a genuine one, and not an artifact 
due to limitation imposed by numerical overflow at low 
temperatures. 

The reentrance found here implies that the long range 
ferromagnetic phase has higher entropy than the low- 
temperature paramagnetic phase. How can we under- 
stand this? A first hint can be obtained by considering 
the behavior of the flow of J(l) close to the upper and 
lower (reentrant) F — P boundary (see Fig. 2). We see 
that J(l) approaches J(l — > oo) — » 0+ monotonously as 
I — * oo close to the upper P— F boundary (curve A). How- 
ever, J(l) swings negative for intermediate length scale 
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(curve B) for all temperatures below the lower (reentrant) 
P— F boundary before eventually approaching the trivial 
paramagnetic fixed point J(l — > oo) = 0. In other words, 
the system establishes short range antiferromagnetic cor- 
relations in the reentrant portion of the phase diagram 
for T < 0.60 and 0.755 < x < 0.855. 
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FIG. 2. Iteration number, Z, dependence of the average 
coupling, J (I) slightly in the paramagnetic phase close to the 
the upper (curve A, x — 0.80, T = 0.70) and lower reentrant 
(curve B, x = 0.80, T = 0.30) P-F phase boundary. 



Interestingly, for the Potts model, a ground state with 
antiferromagnetic correlations in presence of random fer- 
romagnetic bonds has lower entropy than a ferromagnetic 
state with random antiferromagnetic bonds. Consider 
three spins, <r\, &i and 03 with ferromagnetic bonds J12 
and J23. If one of the bond is, instead, antiferromagnetic, 
(T2 becomes an idle and entropy-carrying spin with zero 
effective average exchange field at T = from ferromag- 
netically aligned <j\ and 173. However, for o\ and (73 anti- 
ferromagnetically aligned via the other &( d_2 ' bonds, o~2 is 
in a unique (non-idle) state for a ferromagnetic J12 bond 
and an antiferromagnetic J23 bond. Consequently, anti- 
ferromagnetically correlated triplets of spins (<7i, 02, 03) 
carry lower entropy in presence of random ferromagnetic 
bonds than a ferromagnetic state with random antifer- 
romagnetic bonds. The antiferromagnetic state also has 
lower energy, {E = — J), as compared to the ferromag- 
netic configuration (E=0). Naively, in order to minimize 
the free energy, F = E — TS, this observation suggests 
that, upon cooling, local antiferromagnetic correlations 
should become more and more favorable since entropy is 
less important at low temperatures. This makes plau- 
sible that the system, at low temperatures, prefers to 
form ferromagnetic domains that are antiferromagneti- 
cally aligned (e.g. on intermediate length scales, as found 
in Fig. 2) rather than to keep the long-range ferromag- 
netic order established at higher temperatures. However, 
and this is an important point, true long-range antifer- 
romagnetic order cannot occur since it is known that the 
lower critical dimension for antiferromagnetic order on 



the 6 = 3 hierarchical lattice is four (d = 4) |27|. Thus, 
for a certain concentration range of random AF bonds, 
reentrant behavior from a ferromagnetic phase to a para- 
magnetic phase with local antiferromagnetic correlations 
can occur. In d = 4 (Fig. 3), the F — ► P reentrance 
dissappears and gives rise, instead, as expected from 
the previous argument, to an F — > AF transition upon 
cooling, where here "AF" refers to the Berker-Kadanoff 
phase charaterized by a nontrivial attractive fixed point 
at nonzero temperature fl27[ . 
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FIG. 3. Temperature, T, concentration, x, of antiferro- 
magnetic bonds phase diagram for the Q — 3 state Potts 
spin glass on a d — 4, b = 3 hierarchical lattice. The 
antiferromagnetic— "BK" phase refers to the antiferromagnet- 
ically ordered phase characterized by a fixed point at finite 
coupling J (I = 00) = -10.94661 [30]. 



IV. CONCLUSION 

In conclusion, we have shown that the Q — 3 Potts 
spin glass model on two and three dimensional hierarchi- 
cal lattices undergoes a ferromagnetic to paramagnetic 
reentrance upon cooling. This reentrance is due: (1) to 
the combination of antiferromagnetically correlated spins 
at low temperatures in the phase "rich" in ferromagnetic 
bonds carrying less entropy than ferromagnetically corre- 
lated spins and, (2) to the fact that the lower critical di- 
mension for antiferromagnetic order for the Q = 3 Potts 
antifcrromagnet on hierarchical lattices is four. Conse- 
quently, reentrance occurs at T > in two and three 
dimensional such lattices. The results presented here 
demonstrate that there is no fundamental reason forbid- 
ding a thermodynamic reentrant phase transition in all 
randomly frustrated systems. Our results open the pos- 
sibility that reentrance might occur in some yet to be dis- 
covered real randomly frustrated systems with Euclidean 
lattices. We hope that our results will stimulate further 
studies in that direction. 
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